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ABSTRACT 

We estimate a characteristic timescale for star formation in the spiral arms of disk galaxies, going 
from atomic hydrogen (H I) to dust-enshrouded massive stars. Drawing on high-resolution H I data 
from The HI Nearby Galaxy Survey and 24 (im images from the Spitzer Infrared Nearby Galaxies 
Survey we measure the average angular offset between the H I and 24/mi emissivity peaks as a function 
of radius, for a sample of 14 nearby disk galaxies. We model these offsets assuming an instantaneous 
kinematic pattern speed, £l p , and a timescale, iHi^24/xm, for the characteristic time span between 
the dense HI phase and the formation of massive stars that heat the surrounding dust. Fitting for 
fl p and tHi^24^m, we find that the radial dependence of the observed angular offset (of the HI and 
24 /im emission) is consistent with this simple prescription; the resulting corotation radii of the spiral 
patterns are typically i? CO r — 2.7 R s , consistent with independent estimates. The resulting values of 
£hii->24 ,urn f° r the sample are in the range 1-4 Myr. We have explored the possible impact of non- 
circular gas motions on the estimate of tffl^24/im and have found it to be substantially less than 
a factor of 2. This implies that a short timescale for the most intense phase of the ensuing star 
formation in spiral arms, and implies a considerable fraction of molecular clouds exist only for a few 
Myr before forming stars. However, our analysis does not preclude that some molecular clouds persist 
considerably longer. If much of the star formation in spiral arms occurs within this short interval 
tein24fim, then star formation must be inefficient, in order to avoid the short-term depletion of the 
gas reservoir. 

Subject headings: galaxies: evolution - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: 
spiral - stars: formation 



1. INTRODUCTION 

I Roberts! (|1969l hereafter R69) was the first to develop 
the scenario of spiral-arm-driven star formation in galaxy 
disks. In this picture a spiral density wave induces gravi- 
tational compression and shocks in the neutral hydrogen 
gas, which in turn leads to the collapse of (molecular) 
gas clouds that results in star formation. This work al- 
ready pointed out the basic consequences for the relative 
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geometry of the dense cold gas reservoir 2 and the emer- 
gent young stars: when viewed from a reference frame 
that corotates with the density wave, the densest part 
of the atomic hydrogen (HI) lies at the shock (or just 
upstream from it), while the young stars lie downstream 
from the density wave. Using HI and Ha as the trac- 
ers of the cold gas and of the young stars, respectively, 
R69 found a qualitative support in the data available at 
the time. In this picture, the characteristic timescale for 
this sequence of events is reflected in the typical angular 
offset, at a given radius, between tracers of the different 
stages of spiral-arm-driven star formation. 

While this qualitative picture has had continued popu- 
larity, quantitative tests of the importance of spiral den- 

2 This paper pre-dates observational studies of molecular gas in 
galaxies. 
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sity waves as star-formation trigger (|Lin fc Shul Il964f) 
and of the timescales for the ensuing star formation 
have proven complicated. First, it has become increas- 
ingly clear that even in galaxies with grand-design spiral 
arms, about half of the star formation occurs in loca- 
tions outside the spiral arms ([Elmegreen fc Elmegreenl 
1986). Second, stars formed from molecular clouds 
(not directly from Hi) and very young star clusters are 
dust enshrouded at first. Moreover, the actual physical 
mechanism that appears to control the rate and over- 
all location of star formation in galaxies is the grav- 
itational instability of the ga s and existing stars (e.g. 
ILL Mac Low, fc Kle ssen 20 051). Stars only fo rm above 
a critical density (|Martin fc Kennicuttj 200 lh which is 
consistent with the predicted iToomrei ( 19641 ) criterion 
for gr avitational instability as generalized by iRafikovl 
(2001). Although, in galaxies with prominent spi- 
ral structure local gas condensations are governed by 
magneto-rotational instabilities — spiral arms are regions 
of low shear where the transfer of angular momen- 
tum i s carr ied out by magnetic fields (jKim fc Ostrikerl 
2002, 2006). Obtaining high-resolution, sensitive maps 
of all phases in this scenario (HI, molecular gas, dust- 
enshrouded young stars, unobscured young stars) has 
proven technically challenging. 

If the star formation originates from direct collapse of 
gravitationally unstable gas, and if the rotation curve 
and approximate pattern speed of the spiral arms are 
known, the geometric test suggested by R69 provides a 
timescale for the end-to-end (from HI to young stars) 
process of star formation. Of course, there are other 
ways of estimating the timescales that characterize the 
evolutionary sequence of the interstellar medium (ISM), 
based on other physical arguments. However, other lines 
of reasoning have led to quite a wide range of varying 
lifetime estimates as discussed below. 

Offsets between components such as CO and Ha 
emission in the disks of spiral galaxies have in- 
deed been observed (IVogel. Kulkarni. fc Scovillel 



conversion timescale can vary from the edge of a molec- 
ular cloud (r ~ 4 x 10 6 yr, n a ~ 10 3 ) to the cen- 
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IScoville et alj 120011 ). 
(120061) remarked that 
the angular separation between the dust lanes and the 
peaks of Ha emissi on found for nearby spiral galaxies 
(e.g. observed by iRobertsI 119691 : iRotd 1 19751 ) implied 
timescales of the order o f 10 M yr. More recently 
lEgusa. Sofue. fc Nakanishil (|2004f ). using the angular 
offset between CO and Ha in nearby galaxies, derived 
^co^Ha =s 4.8 Myr. 

Observationally, the HI surface density is found to 
correlate well with sites of star formation and emission 
from molecular clouds (|Wong fc Blitzl [2002; Ke nnicuttl 
1998). The conversion timescale of H Ii— >H2 is a key is- 
sue since it determines how well the peaks of H I emission 
can be considered as potential early stages of star for- 
mation. H2 molecules only form on dust grain surfaces 
in dusty regions that shield the molecules from ionizing 
UV photons. Their formation facilitates the s ubsequent 
build ing up of more complex molecules (e.g. IWiiliamsl 
2005). Within shielded clouds the conversion timescale 
H Ii— >H2 is given by th^ ~ 10 9 /np yr, where rip is the 



proto n de nsity in cm (Hollenbach fc Salpeterl 119711: 
Jural[l975l : [Goldsmith fc Lill2005t iGoldsmith. Li. fc Krcol 
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is higher. Local turbulent compression can further en- 
hance the local density, and thus dec rease the conversion 
timescale (jGlover fc Mac Lowll2007l ). Thus, even short 
cloud-formation timescales remain consistent with the 
H Ii— *H2 conversion timescale. 

Th e subsequent evolution (see e.g. iBeuther et al.l 
120071 for a review) involves the formation of cloud 
cores (initially starless) and then star cluster forma- 
tion through accretion onto protostars, which finally be- 
come main sequence stars. High-mass stars evolve more 
rapidly than low-mass stars. Stars with M > 5Mq 
reach the main sequenc e quickly, in less than 1 Myr 
([Hillenbrand et alJ FiMl IPalla fc Stahlen fl999l ). while 
they are still deeply embedded and actively accreting. 
The O and B stars begin to produce an intense UV 
flux that photoionizes the surrounding dusty environ- 
ment within a few Myr, and subsequen tly become op- 
tically visible {rE ronson fc Telescq|1986|). 

A different scenario is sugg ested by ([2002D . 111 

which young stars in the disks of galaxies produce HI 
from their parent H2 clouds by photodissociation. Ac- 
cording to this scenario, the HI should not be seen fur- 
thest upstream in the spir al arm, but rather b etween the 
CO and UV/Ha regions. I Allen et all (jl986[ ) indeed re- 
port observation of H I downstream of dust lanes in M83. 

Several lines of reasoning, however, point toward longer 
star-formation timescales an d molecular cloud lif e times , 
much greater than 10 Myr. IKrumholz fc McKed (|2005l ) 
conclude that the star-formation rate in the solar neigh- 
borhood is low. In fact, they point out that the star- 
formation rate in the solar neighborhood is ~ 100 times 
smaller than the ratio of the masses of nearby molec- 
ular clouds to their free-fall time Mmc/ t s, which also 
indicates the rate of compression of molecular clouds. 
Individual dense molecular clouds have been argued to 
stay in a fully mol ecular state for about 10-15 Myr be- 
fore their collapse (Ta ssis fc Mouschovias! 120041 ) . and to 
transform about 30% of their mass into stars in > 7 Tg 
(~ 10 6 yr e.g. considering the mass of the Orio n Nebula 
Cluster, ONC I Tan. Krumholz. fc McKed l200l . Large 
molecular clouds have been calculated to survive 20 to 
30 Myr before being destroyed by the stellar feedback 
by IKrumholz. Matzner. fc McKed ([20 061 . Based on ob- 
servations. IPalla fc Stahlerl ([1999L l2000f i argue that the 
star formation rate in the ON C was low 10 7 yr a go, and 
that it increased only recently. iBlitz et all (|2007f ). using 
a statistical comparison of cluster ages in the Large Mag- 
ellanic Cloud (LMC) to the presence of CO, found that 
the lifetime for giant molecular clouds is 20-30 Myr. 

However, other studies conclude t hat the timescales 
for star-formation are rather short. iHartmarml (2003) 
pointed out that the Palla-Stahler model is not con- 
sistent with observations since most of the molecu- 
lar clouds in the ONC are forming stars at the same 
high rate. The stellar age or the age spread in 
young open stellar clusters is not necessarily a use- 
ful constraint on the star-formation timescale: the 
age spread, for example, may result from indepen- 
dent and n on-sim ult aneous bursts of star-formation 
Elmegreenl I2000D . iBallesteros-Paredes fc Hartmannl 
2007) pointed out that the molecular cloud life- 
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time must be shorter than the value of tmc — 

10 Myr suggested by iMouschovias. Tassis. fc Kund 

(2006). Also subsequent star formation must proceed 
very q uickly, within a few Myr (|Vazquez-Semadeni et al.l 
20051: iHartmann. Ballesteros-Paredes. fc Berginl 1200 If ). 
Prescott et al.1 (|2007f ) found strong association between 
24 /xm sources and optical HII regions in nearby spiral 
galaxies. This provides constraints on the lifetimes of 
star-forming clouds: the break out time of the clouds 
and their parent clouds is less or at most of the same 
order as the lifetime of the HII regions, therefore a few 
Myr. Dust and gas clouds must dissipate on a timescale 
no longer than 5-10 Myr. 

In conclusion, all the previous studies listed aim to es- 
timate the lifetimes of molecular clouds or the timescale 
separation between the compression of neutral gas and 
newly formed stars. Most of these studies are based on 
observations of star-forming regions both in the Milky 
Way and in external galaxies, and in all cases the de- 
rived timescales lie in a range between a few Myr and 
several tens of Myr. 

In this paper, we examine a new method (§ [2]) for es- 
timating the timescale to proceed from HI compression 
to star formation in nearby spiral galaxies. We com- 
pare Spitzer Space Telescope/MIPS 24 /im data from 
the Spitzer Near Inf rared Galaxies Survey (SINGS; 
iKennicutt et alj 120031 ) to 21 cm maps from the H I 
Nearby Galaxy Survey (THINGS; I Walter et al.ll200l . 
The proximity of our targets allows for high spatial 
resolution. In § [3] we give a description of the data. 
The MIPS bands (24, 70 and 160 ^m) are tracers of 
warm dust heated by UV and are therefore good indi- 
cators of recent star-formation activity (see for example 
iDale et aT1l2005[ ) . We used the band with the best resolu- 
tion, 24 yitm, which has been recognized as the best of the 
Spitzer bands for tracing star f ormation (jCalzetti et al.l 
[20051 [20071 [Prescott et al][200l : the 8 /im. Spitzer /IRAC 
band has even higher resolution but is contaminated by 
PAH features that undergo str ong depletion in the pres- 
ence of intense UV radiation (Dwek 2005; Smit het al.l 
l2007i ). In §g] we describe how we use azimuthal cross- 
correlation to compare the H I and 24 /im images and 
derive the angular offset of the spiral pattern. This al- 
gorithmic approach minimizes possible biases introduced 
by subjective assessments. We describe our results in §[5] 
where we derive iffl^24^m for our selection of objects. 
Finally, we discuss the implications of our results in § [5] 
and draw conclusions in § [7] 

2. METHODOLOGY 

The main goal of this paper is to estimate geometri- 
cally the timescales for spiral-arm-driven star formation 
using a simple kinematic model, examining the R69 ar- 
guments in light of state-of-the-art data. Specifically, we 
set out to determine the relative geometry of two tracers 
for different stages of star- formation sequence in a sample 
of nearby galaxies, drawing on the SINGS and THINGS 
data sets (see §[3]): the 24 ^m and the HI emission. 

While the angular offset between these two tracers is 
an empirical model- independent measurement, a conver- 
sion into a star-formation timescale assumes (a) that 
peaks of the H I trace material that is forming molec- 
ular clouds, and (b) that the peaks of the 24 /j,m emis- 
sion trace the very young, still dust-enshrouded star clus- 



ters, where their UV emission is absorbed and re-radiated 
into the mid- to far-infrared wavelength range (~5 fj,ra to 
~500/im). The choice of these particular tracers was mo- 
tivated by the fact that they should tightly bracket the 
conversion process of molecular gas into young massive 
stars, and by the availability of high-quality data from 
the SINGS and THINGS surveys. Note that a num- 
ber of imaging studi es in the near-IR have shown (e.g. 
iRix fc Zari tskv 1995|) that the large majority of luminous 
disk galaxies have a coherent, dynamically relevant spi- 
ral arm density perturbation. Therefore, this overall line 
of reasoning can sensibly be applied to a sample of disk 
galaxies. 

We consider a radius in the galaxy disk where the spiral 
pattern can be described by a kinematic pattern speed, 
fl p , and the local circular velocity v c (r) = il(r) x r. 
Then two events separated by a time £hi>->24 ^mwill have 
a phase offset of 

A(f>{r)^(n(r)-n p )t m ^ 24 ^ m , (1) 

where £hih->24 denotes the time difference between two 
particular phases that we will study here. If the spiral 
pattern of a galaxy indeed has a characteristic kinematic 
pattern speed, the angular offset between any set of trac- 
ers is expected to vary as a function of radius in a char- 
acteristic way. Considering the chronological sequence, 
defining the angular phase difference A<p = 024 ~ 0hi 
and adopting the convention that (f> increases in the di- 
rection of rotation, we expect the qualitative radial de- 
pendence plotted in Figure [1] A<fi > where the galaxy 
rotates faster than the pattern speed, otherwise A(j> < 0. 
Where fl(R COI ) = f2 p , at the so-called corotation radius, 
we expect the sign of A<fi to change. 

In practice, the gaseous and stellar distribution is much 
more complex than in the qualitative example of Fig- 
ure [TJ since the whole spiral network, even for galaxies 
where the spiral arms are well defined such as in grand- 
design galaxies, typically exhibits a full wealth of smaller 
scale sub-structures both in the arms and in the inter- 
arm regions. The optimal method to measure the angu- 
lar offset between the two observed patterns is therefore 
through cross-correlation (§ [3J. We treat the timescale 
imi- > 24 (im and the present-day pattern speed fl p as global 
constants for each galaxy, although these two parameters 
might, in principle, vary as function of galactocentric 
radius. Note that we need not to rely on the assump- 
tion that the spiral structure is quasi-stationary over ex- 
tended periods, t > idyn- Even if spiral arms are quite 
dynamic, continuously forming and breaking, and with a 
pattern speed varying with radius, our analysis will hold 
approximately. 

3. DATA 

The present analysis is based on the 21 cm emission 
line maps, a tracer of the neutral atomic gas for the 
14 disk galaxies listed in Table [TJ which are taken from 
THINGS. These high-quality NRAO 3 Very Large Array 
observations provide data cubes with an angular resolu- 
tion of ~ 6" and spectral resolution of 2.6 or 5.2 km s _1 . 
Since the target galaxies are nearby, at distances of 3- 
10 Mpc, the linear resolution of the maps corresponds 

3 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 1. — Schematic geometry adopted to derive the azimuthal 
phase difference (<j>24 — <t>m){r) = A0(r) between the Hi and 
the 24 (im emission, with rfr increasing in the direction of rotation. 
The sketch shows part of a face-on galaxy rotating anti-clockwisc, 
with the center as indicated. The solid curved lines represent the 
two components within one spiral arm, namely the H I and the 
heated dust. The angular separation between the two components 
is exaggerated for clarity. We measured the deprojected phase 
difference A0(r) at a given radius. Inside corotation, -R CO r, the 
material is rotating faster than the pattern speed and the 24 fim 
emission lies ahead of the Hi (</>hi < <^24um)' At corotation the 
two patterns coincide, and outside R CO r the picture is reversed since 
the pattern speed exceeds the rotation of the galaxy. 

to 100-300 pc. The H I data cubes of our target galax- 
ies are complemented with near-IR images, which are 
public data. In particular, the majority of the THINGS 
galaxies (including all those in Table QJ have also been 
observed within the framework of the SINGS and we 
make an extensive use of the 24 jxm MIPS images (see 
§ 14. 2p . Figure [5] illustrates our data for two of the sam- 
ple galaxies, NGC 5194 and NGC 2841. The 24 //m 
band image is shown in color scale, and the contours 
show the HI emission map. To obtain the exponential 
scale length of the stellar disk (see § SJ), we use 3.6 /zm 
Infrared Array Camera (IRAC) images when available, 
otherwise we use H band imag es taken from the Two 
Micron All Sky Survey (2MASS: [Jarrett et al.|[2003l) . To 
check the consistency of our results, we use CO maps 
from the Berkeley-Illinois- Marylan d Aaaociatiqn Sur vey 
Of Nearby Galaxies f BIM A- S ONG IHelfer et ai][2003l) for 
some of our target galaxies. 

4. ANALYSIS 

All analysis in this paper started from fully reduced im- 
ages and data cubes. On this data we carry out two main 
steps. First, we derive the rotation curve v c (r) of the HI 
and the geometrical projection parameters of the galaxy 
disk, and use these parameters to deproject the maps of 
the galaxies to face-on orientation (see Table [T]). Second, 
we sample the face-on maps in concentric annuli. For 
each annulus we cross-correlate the corresponding pair 
of HI and 24 ^m fluxes, in order to derive the angu- 
lar offset between the HI and the 24 /im patterns as a 
function of radius. 

For three of the galaxies listed in Table □ (NGC 628, 
NGC 5194, and NGC 3627), we also measure the angular 
offset between the CO and 24 ^m emission maps. If the 
ISM evolves sequentially from atomic into molecular gas, 
and then subsequently initiates the formation of stars, 
considering the kinematics expressed in Eq. [T] we expect 



the CO emission to lie in between the HI and the 24 /*m. 

4.1. Analysis of the HI Kinematics 

For each object we apply the same general approach: 
we first perform adaptive binning of the HI data cube 
regions with low sign al-to-noise (S/N) ratio us ing the 
method described by ICappellari fc Copinl (|2003h . From 
the resulting spatially binned data cubes we fit the 21 cm 
emission lines with a single Gaussian profile and use the 
parameterization to derive (1) the line-of-sight velocity 
map v(x,y), given by the line centroid, and (2) the flux 
maps fj,o(x,y) = a(x,y)/(^/27r a(x,y)), where a and a 
are the Gaussian peak amplitude and width, respectively. 
Since we do not need to derive the rotation curve with 
high accuracy for the purpose of this paper, we limit our 
model to a co-planar rotating disk with circular orbits 
described by 

v(x,y) — v sys + v c (r) sinicos?/>, (2) 

where v(x,y ) is the observed velocity map along the line 
of sight (see iBegemanl [1989) . For simplicity, we assume 
here that the orbits are circular, though we address the 
issue of non-circular motions in § 15.41 By x 2 minimiza- 
tion fitting 4 of the model function in Eq. [5] to the ob- 
served velocity map v(x,y), we obtain the systemic ve- 
locity v S ys = const, the inclination i and the position 
angle (PA.) of the geometric projection of the disk on 
to the sky. Here, ip is the azimuthal angle on the plane 
of the inclined disk (not the sky) and is a function of i 
and P.A. The line where "0 = denotes the orientation 
of the line of nodes on the receding side of the disk. The 
kinematic center (xo,ya) is fixed a priori and is defined 
as the central peak of either the IRAC 3.6 /nm or the 
2MASS H band image. The positions of the dynami- 
cal centers used her e are c onsistent with those derived in 
iTrachternach et all (|2008h . We parameterize the depro- 
jected rotation curve v n with a f our-parameter arctan-likc 
function (e.g., iRix et al.|[l997l ) 

v c (r) = v (1 + xf (1 + x-*)- 1 '-* , (3) 

where x — r/rQ. Here, ro is the turn-over radius, vq is the 
scale velocity, 7 determines the sharpness of the turnover 
and (3 is the asymptotic slope at larger radii. 

The values for the projection parameters i, and P.A., 
the systemic velocity v sys , and the asymptotic velocity 
that have been obtained applying the approach described 
above, are consistent with the values reported in Table [TJ 
From the maximum value of Eq. [3] we obtain the maxi- 
mum rotational velocity v m&Xl which is listed for all the 
sample galaxies in Table [TJ 

4.2. Azimuthal Cross- Correlation 

The central analysis step is to calculate by what angle 
A<p the patterns of HI and 24 /im need to be rotated 
with respect to each other in order to best match. We 
use the kinematically determined orientation parameters, 
i and P.A., to deproject both the HI and 24 /jm images 
to face-on. To estimate the angular offset A(f> between 
the two flux images at each radius, we divide the face-on 

4 The fitting has been performed with the mpfit IDL routine 

found at the URL: 

http://cow.physics.wisc.edu/~craigm/idl/fitting.html 
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TABLE 1 

THINGS and SINGS Target Galaxies 



Note. 



Obj. Name 


Alt. Name 


R25 (') 


Rs (') 


Band 


i (°) 


P.A. (°) 


D (Mpc) i) max (his- 1 ) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGO 2403 




9.98 


1.30* 


H 


63 


124 


3.22 


128 


NGC 2841 




3.88 


0.92* 


3.6 


74 


153 


14.1 


331 


NGC 3031 


M81 


10.94 


3.63 ±0.2 


3.6 


59 


330 


3.63 


256 


NGC 3184 




3.62 


0.92 ±0.09 


H 


16 


179 


11.1 


260 


NGC 3351 




3.54 


0.86 ±0.03 


3.6 


11 


192 


9.33 


210 


NGC 3521 




4.8 


0.74 ± 0.02 


3.6 


73 


340 


10.05 


242 


NGC 3621 




5.24 


0.80* 


H 


65 


345 


6.64 


144 


NGC 3627 


M66 


4.46 


0.95* 


3.6 


62 


173 


9.25 


204 


NGC 5055 


M63 


6.01 


1.16 ±0.05 


H 


59 


102 


7.82 


209 


NGC 5194 


M51 


3.88 


1.39 ±0.11 


H 


42 


172 


7.77 


242 


NGC 628 


M74 


4.77 


1.10 ±0.09 


H 


7 


20 


7.3 


220 


NGC 6946 




5.35 


1.73 ±0.07 


H 


32.6 


242 


5.5 


201 


NGC 7793 




5.0 


1.16 ±0.05 


H 


50 


290 


3.82 


109 


NGC 925 




5.23 


1.43* 


3.6 


66 


286 


9.16 


121 


Semi-major 


axis of the 


25 mag 


arcsec 2 isophote 


in the 


B band 


obtained from 


the LEDA 



database (URL: 

|http: / /leda. univ-lyonl.fr/} ; (2) exponential scale length derived in this paper as described in § [4] using either IRAF or galfit (values 
tagged with *), where the error bars are 8R S /R S < 1 % ; (3) image band used (2MASS H or IRAC 3.6 jim band) to derive R s \ (4) and (5) 
kinematic inclination and P.A., respectively; (6) adop ted distance; (7) maximum amplitude of the rotatio n velocity corrected fo r inclination 
obtained from the rotation curve v c derived in £j 14.11 The values in Columns (4) to (6) are adopted from dc Blok ct al. (2008). 





Fig. 2. — The 24 /im band image is plotted in color scale for the galaxies NGC 5194 (left) and NGC 2841 (right); the respective Hi 
emission map is overlayed with green contours. 

Expanding the argument of the sum in Eq. [4] one ob- 
tains that x 2 W is independent of the terms J2k x 1 an< ^ 
Y), Uh-i, and x 2 is minimized by the maximization of 

cc x , y (£) = ^ \. x k Vk-l] ' ( 6 ) 

k 

which is defined as the CC coefficient. Here we used the 
normalized CC 

Efc ii x k ~ x) (y k -i - y)] 



images into concentric rings of width ^5" and extract 
the flux within this annulus as a function of azimuth. 
We then use a straightforward cross-correlation (CC) to 
search for phase lags in fm(4>\ r ) versus f24 iim(4>\ r ) ■ I n 
general, the best match between two discrete vectors x 
and y is realized by minimizing as a function of the phase 
shift I (also defined as lag) the quantity 



X 2 

a x, y 



E 

k 



[xk - vk-tY 



(4) 



(7) 



where the sum is calculated over all the N elements of x 
and y with k = 0, 1, 2, N — 1. Specifically here, for a 
given radius r = f we consider for all discrete values of 
azimuth <f>: 

Xk = fm((f>k\f) and yk-i = hinm((j>k-i\r)- (5) 



vEfeO^fe ~ x ) 2 T.k(Vk - y) 2 
where x and y are the mean values of x and y respec- 
tively. Here the slow, direct definition has been used and 
not the fast Fourier transform method. The vectors are 
wrapped around to ensure the completeness of the com- 
parison. With this definition the CC coefficient would 
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have a maximum value of unity for identical patterns, 
while for highly dissimilar patterns it would be much 
less than 1. We apply the definition in Eq. [7] using the 
substitutions of Eq. [5] to compute the azimuthal CC co- 
efficient cc(£) of the HI and the 24 /im images. The best 
match between the H I and the 24 /zm signals is realized 
at a value £ max such that cc(^ max ) has its peak value. 
Since the expected offsets are small (only a few degrees) 
we search the local maximum around £ ~ 0. The method 
is illustrated in Figure [31 which shows that cc(£) has sev- 
eral peaks, as expected due to the self-similarity of the 
spiral pattern. 

We consider a range that encompasses the maximum 
of the cc(£) profile, i.e. the central ~100-150 data points 
around £ max . This number, depending on the angular 
size of the ring, is dictated by the azimuthal spread of 
the spiral arms and the number of substructures (e.g., 
dense gas clouds, star clusters, etc.) per unit area. This 
corresponds, for example, to a width in £ of a few tens of 
degrees at small radii (~ 1'), depending on the distance 
of the object, and a range in width of t decreasing linearly 
with the radius. We interpolate cc(£) around £ max with 
a fourth-degree polynomial using the following approx- 
imation: cc(£) ~ Pi(£) = J2n=o an ^ n an< l calculate nu- 
merically (using the Python 5 package scipy . optimize) 
the peak value at £ m&K , P4(^ m ax)- By repeating the pro- 
cedure for all radii, the angular offset HL-^24 (xm results 
in A<j>(r) — — £ max (r). The direction or equivalently 
the sign of the lag £ max between two generic vectors x 
and y depends on the order of x and y in the defini- 
tion of the CC coefficient. Note that cc XtV (£) in Eq. [6] 
is not commutative for interchange of x and y, being 



the two patterns best 



match at zero azimuthal phase shift. The error bars for 
5£ max (r) have been evaluated through a Monte Carlo ap- 
proach, adding normally distributed noise and assuming 
the expectation values of £ max and <W max as the mean 
value and the standard deviation, respectively, after re- 
peating the determination N = 100 times. 

Our analysis is limited to the radial range between 
low S/N regions at the galaxy centers and their outer 
edges. In the HI emission maps the S/N is low near the 
galaxy center, where the H I is converted to molecular H2, 
whereas for the 24 /im band the emission map has low 
S/N near R25 (and in most cases already at ~ 0.8 i?2s)- 
Regions with S /N < 3 in either the H I or 24 /im images 
have been clipped. We also ignore those points £ max with 
a coefficient cc(£ max ) lower than a threshold cc ~ 0.2. We 
further neglect any azimuthal ring containing less than a 
few hundred points, which occurs near the image center 
and near i?25- The resulting values A<j)(r) are shown in 
Figure 0J 

4.3. Disk Exponential Scale Length 

We also determine the disk exponential scal e length R s 
for ou r sample using the galf it 6 algorithm (jPeng et all 
120021) . In particular, we fit an exponential disk profile 
and a de Vaucouleurs profile to either the IRAC 3.6 /im 
or to the 2MASS H band image. As galf it underesti- 
mates the error on R s (as recognized by the author of the 

5 http://www.python.org 

6 Found at URL: 
http://zwicky.as.arizona.edu/~cyp/work/galfit/galfit.html 



algorithm), typically 5R S /R S < 1%, we ther efore also use 
the IRAF task ellipse (|Jedrzeiew ski 1983) to derive the 
radial surface brightness profile and fit R s . After testing 
the procedure on a few objects, we note only small dif- 
ferences (of the order of the error bars in Table [J) when 
deriving R s from the H band and the 3.6 /im band. 

5. RESULTS 

5.1. Angular Offset 

With the angular offset A</>(r) = (</f>24p m — ^>m)(r), 
where <f> increases in the direction of rotation, and the 
rotation curve v c (r) for each radial bin, we can rewrite 
Eq.[T]as 

'v c (r) 



A(f>(r) 



• 24 |imi 



(8) 



where f2(r) = v c /r. Since fi(r) > fl p inside the coro- 
tation radius i? C or, and 0(r) < il p outside corotation, 
we expect A<j)(r) > for r < R cm and A(f>{r) < 
for r > R COI . At corotation, where A</>(i? cor ) = 0, the 
two components H I and 24 /im should have no system- 
atic offset. We assume iHiH^24/im and fi p to be constant 
for any given galaxy, and that all the spirals are trail- 
ing, since the only spiral galaxies known to have a lead- 
ing pattern are NGC 3786, NGC 5426 and NGC 462 2 
(jThomasson et all 119891 : iBvrd. Freeman. & Butal l2002f ). 
By x 2 fitting the model prediction of Eq. [8] to the mea- 
sured angular offsets A<p(r) in all radial bins of a galaxy, 
we derive best-fit values for tni^24 and fl p . 

The A(j){r) data and the resulting best-fits are shown 
in FigureJH with the resulting best-fit values listed in Ta- 
ble[2] In Figure[4]we plot for all objects the radial profile 
of the angular offsets A<j>{r). The solid line represents the 
best fit model proscribed by Eq. [H The square symbols 
in the plot represent the fitted data points from § 14.21 
Looking at the ensemble results in Figure [4] two points 
are noteworthy: (1) the geometric offsets are small, typ- 
ically a few degrees and did need high-resolution maps 
to become detectable, (2) the general radial dependence 
follows overall the simple prescription of Eq.[8]quite well. 

5.2. <Hi^24pm and R COT 

Because A<p(r) is consistent with (and follows) the pre- 
dictions of the simple geometry and kinematics in Eq. [8l 
the procedure adopted here turns out to be an effec- 
tive method to derive the following: (1) the time lag 
£ffli-^24/jm, which should bracket the timescale needed to 
compress the molecular gas, trigger star formation, and 
heat the dust; it therefore represents also an estimate 
for the lifetime of star-forming molecular clouds; (2) the 
kinematic pattern speed f2 p of the galaxy spiral pattern 
and, equivalently, the corotation radius i? CO r- 

We now look at the ensemble properties of the re- 
sulting values for tm M 24/im and i? C or- The scatter of 
the individually fitted A(f> points is significantly larger 
than their error bars (as shown in Fig. [4]), which may 
be due to the galactic dynamics being more complex 
than our simple assumptions. For example, the pat- 
tern speed may not be constant over the entire disk 
or there may be multiple corotation radii and pattern 
speed s, as, for instance, foun d by numerical simula- 
tions dSellwood fc S parkc 19881) and observed in external 
galaxies ( Hernandez et al.ll2005l ). 
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Fig. 3. — Representative examples for the determination of the azimuthal Hl-24 fim offset: shown is the cross-correlation cc(£) of the 
two functions fmiifiky?) an d f2i pim(4>k— 1> as a function of azimuth offset <f> at a fixed radius f. The present example shows the cc(£) 
profile calculated for NGC 628 at r ~ 2' (~4.2 kpc) and NGC 5055 at r ~ 2.6' (~6.4 kpc), in the left and right columns, respectively. 
Top panel: the cc(l) profile in the entire range [—180°, 180°]; bottom panel: a zoom of the range [—30°, 30°]. We considered an adequate 
range greater than the width of the cc(£) profile around £ ma x and interpolated cc(l) locally (~ ±10° in the two example plots) with a 
fourth-degree polynomial pi(£) = X] n =o a " ant ^ calculated numerically the peak value £max- The bottom panel shows the fit residuals 
overplotted around the zero level. 



TABLE 2 

Characteristic Timescales im^24Mm and Pattern Speed f2 p Resulting from a \ 2 fit of the Observed Angular Offset via 

Equation Q] 



Obj. Name 


Alt. Name 


*HI>-^24 |im 




RcOY / R-S 


n P 






[Myr] 


[km s — 1 kpc -1 ] 




[km s -1 kpc -1 ] 


NGC 2403 




1.4 ±0.5 


30 ±4 


2.8 ±0.3 




NGC 2841 




4.4 ±0.5 


42 ± 2 


2.8 ±0.1 




NGC 3031 


M81 


0.5 ±0.3+ 


27 ± 13+ 


2.3 ± 1.4 


24 a 


NGC 3184 




1.8 ±0.4 


38 ±5 


2.3 ±0.5 




NGC 3351 


M95 


2.2 ±0.3 


38 ±3 


2.3 ±0.4 




NGC 3521 




2.9 ±0.4 


32 ± 2 


3.6 ±0.2 




NGC 3621 




2.3 ± 1.3+ 


31 ± 11+ 


2.8 ± 1.0 




NGC 3627 


M66 


3.1 ±0.4 


25 ± 4 


3.0 ±0.5 




NGC 5055 


M63 


1.3 ±0.3 


20 ±5 


3.8 ± 1.5 


30 - 40 b 


NGC 5194 


M51 


3.4 ±0.8 


21 ±4 


1.5 ±0.6 


38 ±7°, 40±8 d 


NGC 628 


M74 


1.5 ±0.5 


26 ±3 


2.2 ±0.4 


32 ± 2 C 


NGC 6946 




1.3 ±0.3+ 


36 ±4+ 


1.7±0.4 


39±9 C , 42±6 d 


NGC 7793 




1.2 ±0.5+ 


40 ± 10+ 


2.5 ±0.9 




NGC 925 




5.7 ± 1.6 


11 ± 1.0 


2.1 ±0.2 


7.7 f 



Note. — The error bars are evaluated via a Monte Carlo method. The timescales tHl>->24/jm listed in this table are summarized 
in the histogram in Figure \5\ The corotation radius to exponential scale radius R cor /R s ratios a re su mmarized in Figure \Ej W e list 
for comparison in the last colum n ot her measurements of the pat tern speed: (a) IWestpfah] 119911): (b) IThornlev fe Mundvl 11997D; (c ) 
IZimmer. Rand, fc McGrawl pOOlh ; fd) IHernandez et al.l (I2004T) : (el ISakhibov fe Smirnovl 12004li ; (f ) Elmegreen. Wilcots. fe Pisanol H1998D . 
Addressing to the discussion on pattern speeds in § 16.21 we tag as bad fit those galaxies (indicated with the f symbol) where the fit is not 
reliable. 



Even though a considerable intrinsic scatter character- 
izes A(f>(r) for most of our sample galaxies, and the error 
bars of tHi^24^m listed in Table [2] are typically > 15%, 
the histogram of the characteristic timescales t^i^2A fim 
in Figure [5] shows overall a relatively small spread for 
a sample of 14 galaxies of different Hubble types: the 
timescales tmi-*24^m occupy a range between 1 and 4 
Myr for almost all the objects. 

The solid curve in Figure [51 representing the prescrip- 



tion of Eq. intersects the horizontal axis at the coro- 
tation radius R COT , which can be formally derived by in- 
verting Eq. [8] at Acj) — 0. We report in Table [2] the ratio 
between i? cor and the exponential scale length R s for each 
object and show this result in Figure Comparisons of 
our pattern speed fl v measurements wit h other method- 
ologies (e.g. lTremaine fe Weinberglll984l ) are listed in Ta- 
ble[2] The differences with our results may arise since the 
Tremaine- Weinberg method assumes the continuity con- 
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Fig. 4. — Radial profiles for the angular offset Hli— >24 (im for the entire sample, obtained by sampling face-on Hi and 24 (im maps 
concentric rings and cross-correlating the azimuthal profiles for each radius. The solid line is the best-fit model to the observed data 
points, denoted by squared symbols, which has been obtained by x 2 minimization of Eq. \E\ the solid curve intersects the horizontal axis 
at corotation (defined as A(f> = 0). The solid and dashed vertical lines indicate the 2.7 Rs — Rcor value and error bars, derived by 
IKranz, Slvz, fc Rixl d2003h . 

dition of the tracer, which may break down for the gas as 
it is easily shocked, it changes state, and it is converte d 
into stars ([Hernandez et al.ll2005t iRand fc Walrm|[2 004). 
or it can be obscured by dust (|Gerssen fc Debattistal 
2007). Comparison of the observed non-axisymmetric 
motions with hydrodynamical models based on the ac- 
tual stellar mass distribution (jKranz. Slvz. fc Rhd l2003) 
have found a characteristic value of R cor /R s ~ 2.7 ± 0.4 
for a sample of spirals. We plot this range with dashed 
vertical lines in Figure [H whereas the solid curve in 
each panel denotes the actual fit value of R COI /R s for 
each object. Figure [6] illustrates the interesting fact 
that the corotation values found in the present paper, 
R C or/Rs — 2.7 ± 0.2, agr ee well with the completely in- 
dependent estimates that IKranz. Slvz. fc Rixl (|2003f ) de- 
rived by a different approach and for a different sample. 
Therefore, a generic value for the pattern speed of the 
dominant spiral feature of R cor /R s ~ 2.7 ± 0.4 seems 
robust. 



t 








t(HI 



2 3 



4 5 
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Fig. 5. — Histogram of the timescales iHl^24pm derived for the 
14 sample galaxies listed in Table [T] from the fits in Figure [4] (also 
listed in Tabled- The timescales range between 1 and 4 Myr for 
almost all galaxies. 
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NGC 628 



Fig. 6. — The best fit for the spiral arm corotation radius (in 
units of exponential scale radii) obtained by inverting the best-fit 
model of Eq. [T] at A0(r) = 0, and evaluating R s via a bulge-disk 
decomposition on either H band or 3.6 fim band images. The 
solid and dash ed horizontal lines represen t the R CO i/Rs = 2.7±0.4 
value found by Kranz, Slyz, & Rix (2003). The galaxies in the plot 
are sorted by asymptotic rotation velocity (see Table [1} , giving 
no indication of a correlation between dynamical mass and the 
Rcoi/Ra ratio. 

5.3. Comparison With CO Data 

If the basic picture outlined in the introduction is cor- 
rect, then the molecular gas traced by the CO, as an 
intermediate step in the star-formation sequence, should 
lie in between and have a smaller offset from the H I than 
the 24 jum does, but in the same direction. To check 
this qualitatively, w e retrieve the BIMA-SONG CO maps 
■1 Heller et al.ll2003h for the galaxies NGC 628, NGC 5194, 
and NGC 3627. For comparison, we derive the angular 
offset between the CO emission and the 24 /im, applying 
the same method described above for the HI. The results 
are plotted in Figure [71 The scarcity of data points (e.g., 
NGC 628) and their scatter, which is typically larger 
than the error bars, make estimates of tco^24^m and 
R CO r rather uncertain. Therefore, we simply focus on 
A0co-24(r) versus A0Hi-24(?"), which is shown in Fig- 
ure [Jj This figure shows that the values of A</>co-24 all 
lie closer to zero than the values of A(/>Hi-24- 

Hence this check shows that the peak location of the 
molecular gas is consistent with the evolutionary se- 
quence where the HI represents an earlier phase than 
the CO. In this picture the HI has a larger spatial sep- 
aration with respect to the hot dust emission, except at 
corotation, where the three components are expected to 
coincide. However, it is clear that higher sensitivity CO 
maps are needed to improve this kind of analysis. 

5.4. Analysis of Non- Circular Motions 

So far we have carried out an analysis that is based 
on the assumption of circular motions. We quantify here 
non-circular motions and determine to what extent their 
presence affects the estimate of the timescales iHi^24^m, 
which scale with A0. In the classic picture (e.g., R69), 
the radial velocity of the gas is reversed around the spi- 
ral shock, so that the material is at nearly the same 
galactocentric radius before and after the shock. Gas in 
galaxies with dynamically important spiral arms does not 
move on circular orbits, though. Shocks and streaming 
motions transport gas inwards, and gas orbits undergo 
strong variations of direction. If a continuity equation for 
a particular gas phase applies, it implies that in the rest 
frame of the spiral arm, the change of relative velocity 
perpendicular to the arm v± is proportional to the arm to 
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Fig. 7. — Comparison of the angular offsets obtained for 
Hli— >24 (im (squares) and CO^>24 fim (triangles) for the galaxies 
NGC 628, NGC 5194, and NGC 3627. CO maps are taken from 
the BIMA-SONG survey. The solid vertical line in each panel indi- 
cates the position of corotation as obtained by \ 2 fitting of Eq.[l]for 
the H I. These results are qualitatively consistent with a temporal 
star-formation sequence HIh^CO^>24 ^m. 

pre -arm mass flux rat io (for a recent illustration in M51 
see lShettv et al.ll2007l hereafter S07). For example, for 
an orbit passing through an arm with mass density con- 
trast of 10, vj_ would drop by the same factor, producing 
a net inward deflection of the orbit. We do note that the 
gas continuity equation may not actually be valid, since 
stars may form, or the gas may become molecular or ion- 
ized. Non-circular motions could modify the simplified 
scheme of Fig. [T] If the orbit is inward bound near the 
arm, the path of the material between the HI and the 
24 /im arm components is larger than that previously 
assumed for circular orbits. We reconsider the scheme 
of Fig. [T] for a non-circular orbit in the frame corotat- 
ing with the spiral pattern as illustrated in Fig. [5J Here, 
the material moves not along a line at constant radius r, 
but along a line proceeding from a larger radius r + dr, 
specifically from the point A (see Fig. [8]) on the H I arm 
toward the point B on the 24 /im arm, where the two par- 
allel horizontal lines denote the galactocentric distances 
r + dr and r. The material departs from A with an angle 
a inwards (if the material were to proceed instead from 
an inner radius, then a is directed outwards). If a is 
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Fig. 8. — Accounting for non-circular motions in the frame coro- 
tating with the spiral pattern, the two horizontal lines denoted 
with r and r + dr represent two galactoccntric distances. Gas (and 
the resulting young stars) move along the line AB oc A<p, proceed- 
ing from the mean location of dense H I to the mean 24 (im peak, 
both denoted with solid lines, where /3 is the H I pitch angle. The 
gas velocity vector is at an angle a with respect to the circle at 
radius r+dr. The line A'B denotes the (biased) measure of A<f>' at 
constant radius r, which we actually estimate in our data analysis, 
and which is shorter than the actual shift value AB oc A0 by a 
factor of cos/3/ cos(/3 + o), as explained in the text. 



large, the measurement of the spatial shift A'B ~ r A<f>' 
in the simple scheme of circular (deprojected) rings no 
longer represents the actual value of the spatial offset 
AB oc Ac/>, but rather is only a lower limit. Consequently, 
the measurement of the timescale 4hii->24 would be un- 
derestimated, since iHi^24 oc Ac/> from Eq. [8] From 
the geometry of the triangle ABA' in Fig. [8] and since 
A'B oc Ac/>' and AB oc AcA, it can be straightforwardly 
shown that 



Acfi(r) = A</>'(r) 



cos (3 



cos(/3 + a) 



where 



tana = 



-Vr 



(9) 



(10) 



r+dr 



v R.( r ) and v'Ar) = v$ — Q p r are the radial and tangen- 
tial velocity components of the gas in the arm frame, 
respectively, and /3 is the H I arm pitch angle defined by 



tan j3 = 



dr 



(11) 



with (3 — > 90° for a tightly wound spiral. Note that the 
radial and the azimuthal dependences of vr and v'i are 
anchored to each other at the spiral arms through Eq.[TTl 
It follows that if the radial and tangential components of 
the gas velocity, and the pitch angle of the H I spiral arms 
are known, the actual value of Acf> can be calculated by 
applying the correction factor 



k(r) 



cos/3 



cos(/3 + a) 



(12) 



to the directly measured quantity Ac/)'. Note from Fig. [5] 
that if a — > 90 — (3, then k(r) — ► oo, but in this case 
also Ac/)' — > 0, resulting in a finite value for Ac/> through 
Eq. [9] and for the timescale, since £fflH^24^m oc Ac6. 



5.4.1. Methodology 

In these circumstances, if spiral arms in galaxies are 
dynamically affecting the gas kinematics, then it is sen- 
sible to test the effects of non-circular motions especially 
in strong arm spiral galaxies from our sample. We select, 
NGC 5194 (M51), NGC 628, and NGC 6946. There- 
fore, we examine the extreme case, calculating along a 
narrow region along spiral arms how much the stream- 
ing motions affect our angular offset measurements, thus 
the derivation of the timescales im^24/im- Subsequently, 
we generalize the procedure accordingly with our specific 
methodology by averaging these effects over the product 
of HI x 24 /zm fluxes, which is the weighting function of 
the cross-correlation. 

The radial and tangential velocity components are diffi- 
cult to separate unambiguously. The radial component is 
accurately measured along the minor axis, where however 
information on the azimuthal component is lost. The 
opposite occurs along the major axis. Therefore, we fol- 
low the prescriptions for non-circular streaming motions 
analysis from § 3 of S07, to separate the velocity compo- 
nents vr and from the observed HI velocity field fobs, 
adopting 

fobs(-R, <t>) = v sys + (vr(R, 4>) sine* + Vj,(R, <j>) coscA) sini 

(13) 

as a generalization of Eq. [2j Knowing vr and v^, then 
we obtain the geometry of the orbits using Equations [10] 
and [TlJ Specifically, the H I column density and velocity 
map are deprojected and resampled into a polar coordi- 
nate system (R, cA) for simplicity, so that the azimuthal 
phase of the spiral arms can be described by a logarith- 
mic spiral ip = (/> arm - 4>q = ln(_R arm /i? ), for a given 
fiducial radius Rq. We extract the observed velocity v t> s 
along logarithmic lines as illustrated in Fig. [9j which best 
represent the arm phase. Though these logarithmic lines 
are drawn by eye, on top of the spiral arms, we obtain 
reasonable results, e.g., for the case of NGC 5194 we find 
a logarithmic slope of 26°, similar to the 21° found by 
S07. Assuming vr and to be constant along equal 
arm phases, we fit Eq. [13] to the observed velocities ex- 
tracted at each arm phase ip S [4>arm, 4>axm + 27r]. The 
fitting gives vr and as a function of ip and therefore 
as a function of radius through Eq. [IT] (see also S07). 
Provided the angles a and from Eq. \W\ and Eq. [TT] we 
straightforwardly calculate the actual value of the angu- 
lar offset Ac/) given Acf>' from circular orbits assumption 

using Eq. 1 and Eq. p] 

Consistently with iGomez fc Coxl (|2002f l and S07, we 
find in the three considered galaxies that the locations 
of the spiral arms coincide with a net drop-off of the 
tangential velocity and negative radial velocity, possibly 
indicating that near the arms the orbits bend inwards. 
We then calculate, using Eq. [12] the correction k(r) at 
the position of the arms, indicated in Fig. [9] where we 
expect the largest variations for vr and v$. Even so, 
we find for the three considered galaxies that k(r) re- 
sults near unity for all radii, except where v'^ approaches 
zero, which does not necessarily coincide with corota- 
tion, but rather where the division in Eq. [TO] diverges 
and the error bars are large. In particular, for NGC 5194, 
where |a| < 20° for all r, the correction k(r) ranges be- 
tween 0.7 and 1.2, and for both NGC 628 and NGC 6946 
k(r) ranges between 0.9 and 1.5, as indicated in Fig. [TOl 
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Fig. 9. — The figure illustrates logarithmic spirals of the form 
if 1 = 0arm — 00 = ln(i? arm / Rq ) where the observed velocity v Q \y a 
is extracted in order to be fitte d to Eq. 1131 and obtain the vr and 
v$ velocity components (see § 15.41 1. The top and bottom pan- 
els represent the projection into a polar coordinates system (R, 0) 
of the observed Hi line-of-sight velocity field and the Hi column 
density map, respectively, for the galaxy NGC 5194. The solid 
lines represent logarithmic spiral arms with phases ij) = a rm and 
= 0arm ± 180°. In the top panel, the gray scale image of the 
velocity field indicates velocity values from —90 km s _1 (dark) 
to +90 km s — 1 (light). The velocity map in figure is not depro- 
jected for inclination effects, which are instead taken into account 
in Eq. 1131 The coordinate = represents the kinematic position 
angle of the galaxy. 

After calculating the corrected offsets A(f> from Eq. [9l 
and fitting Eq. [8] to the values A(j>, we find that the 
timescale iffl>-^24^m and the pattern speed fl p do not 
change significantly — the differences are below the error 
bars for the three galaxies. Note that the data points 
where v'^ ~ 0, that is where Eq. [10] diverges, are not ex- 
cluded from the fit. The results of these fits are listed 
in Table [3] and plotted in Fig. [TO] We also find that the 
radial displacements, 

dr ~ r tan[A<//(Y)] k(r) sin[a(r)], (14) 

are typically as small compared to the radial steps of 
A<f>', i.e. |dr| < 70 pc for all r for NGC 5194. In 
§ 14.21 we calculated A</>' through Eq. [7J hence not only 
within the spiral arms, but also as intensity-weighted 
mean across all the azimuthal values. If we were to calcu- 
late the angle-averaged value A<j>' from the average (vr) 
and {v$) weighted by the product HI x 24 //m, which is 
the weighting function of the cross-correlation, we find 
a ~ for all radii and a correction k(r) even closer to 
unity than over a region limited to the arms. With this 
approach we obtain 0.95 < k(r) < 1.05 for NGC 5194. 

After exploring the two extreme cases, (1) simple 
model of circular orbits and (2) streaming motions near 
spiral arms for galaxies with prominent spiral structure — 
where these effect are supposed to be the largest — we 
find that the implied timescales £m>-^24/jm do not vary 
significantly. By estimating the streaming motions in 
three galaxies from our data set, we find that the cor- 
rection k(r) which we must apply to the angular offset 
measurements in the scheme of circular orbits is gener- 
ally near unity. Non-circular motions do not greatly af- 
fect the offset measurements A<j> for the galaxies with the 
most prominent spiral arms of our data set, where we ex- 
pect indeed the highest deviations from circular orbits, 
suggesting that, since tHL->24/im cx A<f>, the timescales 
iHii-^24/jm will not vary by more than a factor of 1.5. 

However, these conclusions should be viewed cau- 
tiously, since (1) we are assuming that vr and can be 
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Fig. 10. — Correction factor k(r) from Eq. 1121 as applied to the 
angular offset measurement for the galaxies NGC 5194, NGC 628, 
and NGC 6946. Bottom panels: the solid curve represents the 
correction factor k(r) = cos/3/ cos(/3 + a) as a function of radius. 
Top panels: the squares denote the angular offset measurements 
A0' of Fig. [4]calculated assuming circular orbits, the circles denote 
A0 after correction, and the solid curve represents the model fit of 
Eq. [8] to the corrected offset values. 

obtained by fitting Eq.[l3l and (2) we do not account for 
beam effects. While the presence of non-circular motions 
produces apparent radial variations of inclination and 
position angle, our estimates are suggesting that these 
effects do not influence much the determination of the 
timescales. Yet, tidal interactions could ensue physical 
variations of inclination and position angle, thus of the 
actual geometry of the observed velocities. Although we 
note that among the sample galaxies only NGC 3031 
(M81), NGC 3627, NGC 5055, and NGC 5194 are af- 
fected by tidal interaction. Moreover, we do not con- 
sider extra-planar motions due to the implied numerical 
difficulties. For instance, to include a vertical velocity 
component v z into Eq. [TS] would introduce a degeneracy 
while fitting Vr, v r f >1 and v z , and a degeneracy in the 
geometry of the motions, rendering the estimates of the 
timescales uncertain. However, after subtracting a cir- 
cular orbit model from the observed velocity field of the 
galaxy NGC 3184, an inspection of the velocity residuals 
reveals deviations from circular motions of ~5-10 km s _1 
amplitude, and about zero near the spiral arms. Consid- 
ering that NGC 3184 is nearly face-on (i = 16°), then 
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TABLE 3 

Characteristic Timescales *hi^24 and Pattern Speed Q p Resulting from a x 2 Fit to the Angu lar Offsets A(j> After 
Correction for Non-circular Motions Following the Prescriptions in Section 15.41 



Obj. Name 


Alt. Name 


*HI>-.24 urn 


S 2p 






[Myr] 


[km s — 1 kpc -1 ] 


NCC 5194 


M51 


3.3 ±0.6 


20 ±3 


NGC 628 


M74 


1.4 ±0.5 


26 ±3 


NGC 6946 




1.1 ±0.3 


36 ±4 



Note. — The best fits and the corrections are plotted in Fig. 1101 In comparison with the values listed in Table [2] <hi^24 urn an d 
differ by less than their corresponding error bars. 



vertical motions should not exceed 5-10 km s _1 . The 
beam deconvolution, on the other hand, would enlarge 
the uncertainties on separating vr and but we rely 
on a resolution size limit which is far below the typical 
thickness of a spiral arm and allows us to resolve fine 
sub-structures within the arms. We also rely on a large 
number of data points — several hundred to several thou- 
sand depending on the galactocentric radius. 

6. DISCUSSION 

By analyzing the angular offsets between HI and 
24 /im in the context of a simple kinematic model, we 
found short timescales, iffl^24^m, as summarized in the 
histogram of Figure [5] and in Table [2l The implied char- 
acteristic timescales for almost all sample galaxies lie in 
the range 1-4 Myr. This result sets an upper limit to the 
timescale for massive star formation under these circum- 
stances, since we are observing the time lag between two 
phases: a) the atomic gas phase, which subsequently is 
compressed into molecular clouds and forms clusters of 
young, embedded, massive stars, and b) the warm dust 
phase, produced by heating from young stars, whose UV 
radiation is reprocessed by the dust into the mid-IR, as 
observed at 24 fxm. For the few objects where there are 
suitable CO data, we checked that this geometric picture 
also holds for the molecular phase. 

6.1. Timescales Derived from Pattern Offsets 

lEgusa, Sofue. fe Nakanishil (|2004D used the same an- 
gular offset technique to compare CO and Ha emis- 
sion. They report timescales <cowHa — 4.8 Myr for 
the galaxy NGC 4254 (not in our sample). The Ha 
traces a later evolutionary stage than the warm dust 
emission, which can indicate the presence of a young clus- 
ter still enshrouded by dust. Therefore, the Ha and the 
dust emission are expected to be separated by the time 
needed to remove the dusty envelope, th ough they are 
observed to be s patia lly well correlated (IWong fc Blitzl 
120021 iKennicuttj (l998fh iPrescott et all (|2007f) found a 
strong association between 24 Lim sources and optical 
HII regions in SINGS galaxies. Also infrared sources 
located on top of older, UV-bright, cl usters that do 
not h ave Ha emission are rare. Since IPrescott et al.l 
(2007) suggest that the break out time from dust clouds 
is short (~1 Myr), we do not e xpect a strong offset. 
lEgusa. Sofue. fc Nakanishil (|2004l ) derived the angular 
offset by subjective assessment of the separation of the 
intensity peaks. They report that this may be a source 
of systematic errors, since they cannot detect by eye, the 
angular phase differences less than a certain threshold 
and in turn, their results would possibly be an upper 
limit. Given these considerations, the timescales derived 
in this paper are likely consistent with the conclusions of 



lEgusa. Sofue. fc Nakanishil (l200l. 

iRotsI (119751) andlG arcia-Burillo. Guelin. fc Cernicharq 
11993D found time lags of ~10 Myr for M81 and M51, 
respectively. In particular, Rots applied the angular off- 
set method to the dust lanes and Ha. This may be in 
part problematic since the dust absorption in the optical 
bands only traces the presence of dust and it is unre- 
lated to the warm dust emission due to star-formation 
onset. Garcia-Burillo et al. measured the spatial separa- 
tion projected on the sky between CO and Ha and not 
t he azimuthal offset. 

lAllenl (|2002f ) has argued that H I is a photodissociation 
product of UV shining on molecular gas, so it should be 
seen b etween the CO and UV/Ha regions. Allen et all 
(1986) observed HI be tween spiral arm dust lanes and 
HII regions. However, [Elmcgrcen| (|2007f) points out that 
there is no time delay between dust lanes and star forma- 
tion: dust lanes may only represent a heavy visual extinc- 
tion effect and may not be connected to star-formation 
onset. Our finding that CO is situated between HI and 
hot dust (i.e. Fig. [7]) stands in confli c t with the predic- 
tions of the model proposed by I Alien! (|2002f ). 

The evolutionary timescales of the ISM phases, espe- 
cially for star formation, are not well constrained. The 
observational results of the last few decades, arrive at dif- 
ferent and, in some cases, controversial conclusions. The 
discrepancies might in part be attributed to effects of 
limited resolution (see § 16 . 3[) . This suggests that higher 
resolution and sensitivity maps, especially for the CO 
emission, are needed to improve the presented technique 
in the future. 

6.1.1. Photodissociation of H2 

We argue that illumination effects of UV radiation 
shining on molecular and dusty regions cannot pho- 
todissociate molecules in order to produce the observed 
peaks of H I — which correspond to typical surface densi- 
ties of several solar masses per pc 2 . First, the mean free 
path of the UV photons is remarkably short, typically 
~100 pc. The presence of dust, particularly abundant in 
disk galaxies, is the main source of extinction in particu- 
lar within spiral arms, where we observe the peak of dust 
emission. Second, none of the galaxies from our data set 
presents prominent nuclear activity, whose UV flux could 
ionize preferentially the inner surfaces of the molecular 
clouds. Also, we exclude that UV radiation from young 
stellar concentrations can ionize preferentially one side 
of the clouds causing HI and CO emissions to lie offset 
with respect to each other, which we instead interpret as 
due to an evolutionary sequence. In fact, if the light from 
young stars effectively produces H I by photodissociation 
of H2, then the neutral to molecular gas fraction is ex- 
pected to increase with star-formation rate per unit area, 
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Ssfr, as a consequence of the increasing UV radiation 
flux. Yet, the HI to H2 ratio decreases with increasing, 
Esfr , as also shown, for example, in iKennicutt et all 
(2007) for the galaxy M51. Moreover, the HI density 
does not vary much as a function of Ssfr (Kc nnicuttl 
fl998h . 

6.2. Pattern Speeds 

The results of the fits in Figure [4] are consistent with 
the existence of a kinematic pattern speed for the consid- 
ered galaxies, and are suggesting that the spiral pattern 
must be metastable at least over a few Myr or, strictly 
speaking, quasi-stationary. If the apparent spiral struc- 
ture seen in young stars were pro duced by stochastic self- 
propagating star form ation fe.g. iGerola fc Seidenl 119781 : 
ISeiden fc Gerolalll982f ) and shear, without an underlying 
coherent mass perturbation, then presumably we would 
not observe the systematic radial variation of the off- 
sets, seen in Figure 0] in particular not that the offset 
changes sign at R ~ 2 — 3i? cxp . Our results, however, 
do not exclude the stochastic star-formation mechanism 
to occur, rather that this is not the dominant trigger of 
star formation. The spiral pattern might be the mani- 
festation of full wealth of modes of propagating density 
waves, which are continuously forming and dissolving gas 
clouds and structures, where, i.e. the azimuthal modes 
for a grand-design spiral are dominant at low orders (e.g., 
m = 2,3,4). If the spiral structure is quasi-stationary, 
then it can be characterized by an instantaneous pat- 
tern speed to first order — at least over a timescale much 
shorter than the orbital time. This is opposed to den- 
sity waves dynamically driven by bars or interaction with 
companions, which could last a few orbital times. Note, 
however, that our analysis holds approximately in both 
possible cases. 

Some of the fits in Figure 2] do not accurately mirror 
the trend of the observed data points, where the scat- 
ter is so large that it renders the interpretation prob- 
lematic. In particular, we could designate as a bad 
fit the results for the galaxies NGC 3031, NGC 3621, 
NGC 7793, and NGC 6946 (Tabled]). We also note that 
the pattern speed obtained in this paper for the galax- 
ies NGC 5055 and NGC 5194 disagree with previous in- 
dependent measure ments (e.g. lThornlev fc MundvllT997l ; 
IZimmer et aLll2004[ ). The fit for NGC 5055 presents large 
error bars due to the large scatter in the azimuthal off- 
sets (see Figure S]), which may explain the differences. 
For NGC 5194 the large difference between our result and 
the Tremaine- Weinberg method prediction could be due, 
as mentioned in § 15.21 to the assumptio n of continuity 
for the gas. In fact, IZimmer et~aT] (|200l find a pattern 
speed f2p ~ 38 km s _1 kpc -1 that is much faster than the 
patt ern speed predicted by the hydrodynamical models 
from lKranz et alj (|2003[ ). fi^ ~ 12 km s" 1 kpc" 1 , using 
i? cor ~ 2.7 R s and R s ~ 1.4'. Moreover, the presence 
of large variations of the azimuthal offsets as a function 
of radius with respect to the smoother fitted curves of 
Figure |4] suggests that the pattern speed may not be 
constant over the entire disk. Instead, the spiral pat- 
tern could be described by more than one pattern speed, 
implying that the delay time is not exactly the same in 
all parts of an individual galaxy, and not necessarily the 
same in all the considered galaxies. However, the off- 



sets measurements and the implied pattern speeds and 
timescales could be statistically pointing toward regions 
of high H I and 24 /im fluxes when using the weighting 
of Equation [3 therefore toward high Esfr, ensuing the 
timescales to be comparable in all cases. 

6.3. Is Star Formation Triggered by Spiral Waves? 

Addressing to our results (§[5] and Fig. [4j) , a further as- 
pect emerging from our analysis concerns the same gen- 
eral behavior displayed by an heterogeneous sample of 
galaxies going from grand-design (e.g. NGC 5194) to 
flocculent morphologies (e.g. NGC 2841 and NGC 5055, 
see also Fig. [3] for an example comparison). Visual ex- 
amination of the IR band images for our sample galaxies 
(e.g. at 3.6 /im) indicates that almost all have two-arm or 
multi-arm coherent spiral arms; for none of these galax- 
ies the spiral structure is so chaotic as to be character- 
ized as flocculent. Moreover, both types of galaxies from 
our sample display comparable integrated star-formation 
rates. Grand-design density waves are not likely to be 
the primary trigger of star formation, since a substantial 
portion of stars are also formed in the inter-arm regions 
(Elme green fc E lmcgrccnl Il986| ) . The concentration of 
young stellar populations near prominent spir al arms is 
rather an effect of kinematics (iRobertsl [l969f). Grand- 
design and flocculent galaxies exhibit the same intrinsic 
self-s imilar geometry (jElmegrcen. Elmeg reen. fc Leitnerl 
l2003h . The difference between these two types of galaxies 
is only dictated by a different distribution of azimuthal 
modes. Large scales structures (low-order modes) are 
the dominant features in grand design. However, the 
short timescales suggest that the physical scales where 
star formation occurs are rather small, since this time 
delay cannot exceed the free-fall time. Weak density 
waves, those described by high-order azimuthal modes, 
not necessarily grand-design modes, are likely to facili- 
tate the growth of super-critical structures which end in 
star-forming events. In galaxies morphologically classi- 
fied as flocculent and grand design, the mechanism that 
triggers star formation must be the same — gravitational 
instability. Direct compression of gas clouds is indeed 
able to locally trigger star formation. 

The sizes of the structures of active star formation, 
however, are not representative of the timescales in- 
volved. Star-forming regions are organized hierarchically 
according to the small-scale turbulent motions of gas 
and stars, where the dyna mical time varies as a func- 
tion of the local scale sizes (lEfremov fc Elmegreenlll998l : 
iBallesteros-Paredes et al.lll999[ ). The improved accuracy 
of recent observational techniques, e.g. THINGS and 
SINGS, allows us to observe smaller and smaller struc- 
tures, which evolve thus more rapidly and are charac- 
terized by shorter timescales. Stellar activity and turbu- 
lence limit the lifetime of molecular clouds by causing the 
destruction of the parent cloud and the cloud dispersal, 
respectively. This lifetime, typically 10 — 20 Myr, drops 
by a factor of ~10 for the star-forming clouds. More- 
over, star formation begins at high rate in only a few 
Myr, appearing in large structures of O-B complexes — 
beads on a string — of few hundred pc scales, and re- 
maining active for ~30-50 Myr, but with a gradually 
decr easing star-formi ng rate, until the complete quench- 
ing (|Elmegreenll2007f ) . The timescales measured here re- 
fer to the very initial phase of star formation, specifically 
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when star-formation rate is the highest — traced by 24 /im 
peak emission. The time delay tm^24^m needs not to 
represent the average time difference between these two 
phases, especially if there is a gradually declining tail 
of star- formation activity, taking longer than tHii->24/iin- 
If much of the star formation were actually occurring 
within our short ^Hi^24^m estimate, then star formation 
would have to be inefficient to avoid conflicts with the 
short-term depletion of the gas reservoir. 

6.4. Can We Rule Out Timescales of 10 Myr? 

Since the 24 /im emission traces the mass-weighted 
star- formation activity, the timescale iffli->24/im mea- 
sures the time for star clusters to form from HI gas, 
but it does not show that all molecular clouds live 
only ^2 Myr. It might show that the bulk of mas- 
sive stars that form in disks has emerged from molec- 
ular clouds that only lived ~2 Myr, though there could 
still be molecular clouds that live an order of magnitude 
longer. The full molecular cloud lif etimes estimated for 
the LMC correspond t o ~10 Myr (jMizuno et all 120011 : 
lYamaguchi et al.ll2001h . while star clusters are formed 
from molecular clouds in only few Myr. Yet, more re- 
cent ob servations suggest much slower evolution. In par- 
ticular, iBlitz et all (|2007l ) propose molecular cloud life- 
times for the LMC as long as 20-30 Myr. In the Milky 
Way this timescale is estimated to be a few M yr (e.g. 
lHartmann. Ballesteros-Paredes. fe B ergin 200^), but the 
examined cloud complexes, such as Taurus and Ophi- 
uchus, have low star-formation rate and masses more 
th an an order of ma gnitude lower than those studied 
by IBlitz et all (|2007t ). Thus, it remains possible that 
our results and these previous arguments are all consis- 
tent. At high density st ar formation also occurs at higher 
rate. IBlitz et all ((20071 ) also find that the timescales for 
the emergence of the first H II regions traced by Ha in 
molecular clouds is inn ~ 7 Myr, which does not exclude 
that tmi > iffl^24/im, but it coul d be problematic wit h 
the short break out suggested bv lPrescott et al.1 (|2007h . 
since the time delay which is required between the onset 
of star formation (as traced by 24 /jm from obscured H II 
regions) and the emergence of Ha emission would need 
to be large. 

In conclusion, we note that the timescale tni^24 ^im re- 
sults from a \ 2 fit to all the data, and it is treated as 
a global constant individually for each galaxy, so it is 
not a function of radius. Globally, all the characteristic 
timescales tHii-^jUtn listed in Table[2]are <4 Myr, except 
one single case (NGC 925). The error bars are also rela- 
tively small: <1 Myr for the majority of the cases. These 
results clearly exclude characteristic timescales iHi>-»24 
of the order of ~10 Myr, even for the highest value 
recorded in our data set which is NGC 925. 

6.5. Theoretical Implications 

The short timescale found here between the peak of 
H I emission and the peak of emission from young, dust- 
enshrouded stars has implications for two related theoret- 
ical controversies. First is the question of whether molec- 
ular clouds are short-lived, dynamically evolving object s 
( Ballester os-Paredes. Hartmann. fe Vazauez-SemadenH 



12001 o r quasi-st a tic ob j e cts evolving over many free-fall 
times (jMatznerl 120021 : IKrumholz. Matzner. fe McKed 
2006). The second, related question is what the rate- 
limiting step for star formation in galaxies is: formation 
of gravitationally unstable regions in the HI that 
can collapse into molecular clouds (lElmegreenl 120021 : 
lKravtsovll200l iLi et al J 120051 120061 lElmegreenll2007lh or 
formation of dense, gravitatio nally unstable cores within 
quasi-stable molecular clouds (IKrumholz fc McKe dl200 



Krumholz. Matzner. fc McKed 120061 : IKrumholz fc Ta 
20071) 



19991 lHartmann et alj 120011 : lElmegreenl 120001: 
Ballesteros-Paredes fc Hartmann! 120071 : lElmegreenl 



The short timescales found here for the bulk of mas- 
sive star formation in regions of strong gravitational in- 
stability appears to support the concept that molecular 
cloud evolution occurs on a dynamical time once gravi- 
tational instability has set in, and that the rate-limiting 
step for star formation is the assembly of HI gas into 
gravitationally unstable configurations. Our work does 
not, however, address the total lifetime of molecular gas 
in these regions, as we only report the separation be- 
tween the peaks of the emission distributions. Molecular 
clouds may well undergo an initial burst of star forma- 
tion that then disperses fragments of molecular gas that 
continues star fo rmation at low ef ficiency for substantial 
additional time (|Elmegreenl 120071 ) . Averaging over the 
efficient and inefficient phases of their evolution might 
give the overall low average values observed in galaxies 
(e.g. IKrumholz fc~TaI3 [2007Tl . 

7. CONCLUSIONS 

We have derived characteristic star-formation 
timescales for a set of nearby spiral galaxies, using 
a simpl e geom etric approach based on the classic 
iRobertd (|1969f ) picture that star formation occurs just 
downstream from the spiral pattern, where gas clouds 
have been assembled into super-critical configurations. 
This derived timescale, tm^i^m, refers to the pro- 
cesses from the densest HI, to the molecular phase, to 
enshrouded hot stars heating the dust. The analysis is 
based on high-resolution 21 cm maps from THINGS, 
which we combined with 24 ^m maps from SINGS. We 
assume that the observed spiral arms have a pattern 
speed ftp. Given the rotation curve, v c (r) — rO(r), this 
allows us to translate angular offsets at different radii 
between the HI flux peaks and the 24 /im flux peaks in 
terms of a characteristic time difference. 

At each individual point along the spiral arm we found 
considerable scatter between the H I and 24 /im emis- 
sion peaks. However, for each galaxy we could ar- 
rive at a global fit, using a cross-correlation technique, 
and derive two characteristic parameters, im>->24/jm and 
i? cor . For our 14 objects we found the general relation 
i? cor = (2.7 ± 0.2) R Sl which is consiste nt with previ- 
ous studies (e.g. lKranz. Slvz. fc R ix 2003) and, more im- 
portantly, we found fmi->24 J ura to range between 1 and 
4 Myr. Even when accounting for uncertainties, at the 
highest peak of star-formation rate timescales as long as 
iHi>->24/am ~ 10 Myr, which have been inferred from other 
approaches, do not appear consistent with our findings. 
At least for the case of nearby spiral galaxies, our analy- 
sis sets an upper limit to the time needed to form massive 
stars (responsible for heating the dust) by compressing 
the (atomic) gas. Therefore, it points to a rapid proces- 
sion of star formation through the molecular-cloud phase 
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in spiral galaxies. If star formation really is as rapid as 
our estimate of im^^m suggests, it must be relatively 
inefficient to avoid the short-term depletion of gas reser- 
voirs. 
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